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Abstract 

Invasive electroencephalograph (EEG) recordings of ten patients suffering 
from focal epilepsy were analyzed using the method of renormalized entropy. 
Introduced as a complexity measure for the different regimes of a dynamical 
system, the feature was tested here for its spatio-temporal behavior in epilep- 
tic seizures. In all patients a decrease of renormalized entropy within the 
ictal phase of seizure was found. Furthermore, the strength of this decrease 
is monotonically related to the distance of the recording location to the fo- 
cus. The results suggest that the method of renormalized entropy is a useful 
procedure for clinical applications like seizure detection and localization of 
epileptic foci. 
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I. INTRODUCTION 



Focal epilepsies are characterized by seizures which originate from a distinct region of 
the brain. The identification of this so called epileptogenic focus is a prerequisite for surgical 
treatment. Thus much attention is paid to the characterization of invasive electroencephalo- 
graph (EEG) recordings of patients suffering from this disease. Competitive techniques like 
high-resolution positron-emission tomography or MRI using specific ligands or metabolites 
are used as alternative to localize the epileptogenic focus in order to avoid invasive EEG 
recordings. For the time being no superiority of either method has been established yet 
Any method will be preferred which has a better spatial and timewise resolution than 
the other with acceptable validity. Within this scenario a new approach in EEC-analysis is 
described in this paper. 

Since there is only little doubt about the nonlinearity of the dynamical system underlying the 
observed time series a broad range of nonlinear analysis techniques has been applied to these 
data. As one of the first results of nonlinear EEC analysis indication of low-dimensional 
chaos was reported by Babloyanz and Destexhe who claimed a noninteger correlation di- 
mension in the EEC recording of a petit mal seizure. Furthermore, Frank et al. obtained 
positive Lyapunov exponents by analyzing two seizure recordings. Similar effects, giving 
indication of chaotic behavior in EEC data, were also described by Freeman and Skarda 
0, Dvorak and Siska 0, Basar and Bullock and Pijn et al. [ p!0[ 



Windowed estimates 



of these measures for piecewise classification of EEC were used by Lehnertz and Elger [IT 
Tirsch et al. [|1^, Pritchard et al. [jl3[ and lasemidis and Sackellares [|14|. Lehnertz et al. pT[ 
analyzed the EEC recordings of 20 patients suffering form unilateral temporal lobe epilepsy. 
The variability of the correlation dimension, estimated for subsequent segments of the EEC, 
was found to be a good indicator of the lateralisation of the epileptic focus. In a case study 
Lerner |]T5[ found the correlation integral itself to be suitable to detect seizure activity in an 
EEC recording. 

By now it is well appreciated that results, obtained by use of these measures, have to be 
interpreted with care. Several investigators 0-18| have especially shown that the value of 
the correlation dimension is influenced by computational as well as recording parameters. 
Theiler et al. |19| found no evidence for low-dimensional chaos when re-examining EEC data 
of 110 patients. The finite correlation dimensions, they obtained for the same data sets in 
an earlier study ||2^, were found to be caused by an artifact of the autocorrelation in the 
oversampled signals. 

A methodologically different approach is based on the examination of neural spike trains 
which might be observed in the EEC of epileptic patients. In a basic implementation the 
statistical properties like mean, variance or skewness of the interspike interval distribution 
are used to characterize the EEC. Although this method neglects the sequence of these in- 
tervals it is proven to be of practical use for some applications |^ . However, for analysis of 
epileptic seizures advanced methods seemed to be necessary, i.e sequence-sensitive methods 
P^^5[. Using a sequence-sensitive complexity measure Rapp et al. found a decrease of 
nonrandom structure of these sequences in focal epileptic seizures induced in rats by appli- 



cation of penicillin p6 



A qualitative characterization of epileptic seizures is given by Heyden et al. [^]. Their 
results, obtained by analyzing EEC data of patients suffering from mesial temporal lobe 
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epilepsy, indicate that the property of reversibihty of a time series can be used to discrimi- 



nate between seizure and non-seizure activity. Casdagh et al. ||28[ report recurrent activity 
to occur in spatio-temporal patterns related to the location of an epileptic focus. 
In the present paper a procedure is described which seems to be appropriate to classify the 
EEG of epileptic patients in an uniform way. The method of renormalized entropy, proposed 
in for classification of the different states of a djTiamical system, is applied to the EEG 
data of ten patients suffering from temporal lobe epilepsy. The method is tested for its abil- 
ity to assign a given segment of an EEG time series to the corresponding neurophysiological 
state of the epileptic patient, i.e. interictal, ictal or posictal phase, as well as its use for 
localization of an epileptic focus. 



II. RENORMALIZED ENTROPY 

Applying Klimontovich's S-theorem pOl to Fourier spectra of scalar time series, Saparin 



et al. introduced renormalized entropy as a complexity measure for the different regimes 
of a dynamical system. Here each regime is represented by the normalized Fourier spectrum 
Si{uj) of one observable Xi{t) which, in formal analogy to classical statistical physics, is 
viewed at as a distribution function of an open system. 

Given a reference distribution Sr{uj), representing the systems state of equilibrium, the 
relative degree of order of the regime described by Si^uj) is determined by comparing the 
entropies of these two distributions under the additional constraint of same mean energy in 
both states. To this end an effective Hamiltonian 

H,ff = -\nSriuj) (1) 

is introduced and the reference spectrum is renormalized ("heated up") into 

Sr{uj) = C(r,)e~^^ = C{T,)Sr{ujy/^\ (2) 

so that holds: 

J Sriuj)Heff{uj)duj = J Si{uj)Heff{uj)du (3) 

and 

' Sr{oj)du = l. (4) 



Here eq. (3) ensures the equality of mean energies of the two states, eq. (4) the normalization 
of the renormalized spectrum. C(Tj) is a normalization factor depending on Tj. 
Renormalized entropy now is given by: 

AH = J Sr{id)\nSr{id)duj — J Si{uj)\nSi{uj)du . (5) 



Applying this method to the different regimes of the logistic map, Saparin et al. ||2^ found 



renormalized entropy to clearly detect all transitions between the different types of periodic 
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behavior as well as the different types of chaos. Kurths et al. [^] and Voss et al. analyzed 
the heart rate variability of patients after myocardial infarction. They report renormalized 
entropy to be a suitable method for the detection of high risk patients threatened by sudden 
cardiac death. 

Because of eq. (1) the method can only be applied to processes which have a purely positive 
spectrum, i. e. chaotic or stochastic processes. Because spectra of such processes have to 
be estimated the spectrum chosen as reference should be of lower energy than each other 
state of the system to avoid "temperatures" Tj less than 1. "Cooling down" an estimated 
spectrum increases the variance of the estimate and therefore the variance of the estimated 
entropy. 



III. ANALYSIS OF THE EEG 

The EEG data analyzed in this study were recorded using chronically implanted subdural 
and intrahippocampal electrodes measuring the local field potential. The obtained signals 
were passed to a multi-channel amplifier system with band-pass filter settings of 0.53Hz — 
85Hz and were written to a digital storage device with a sampling interval of At = 5.76 ms 
per channel. Fig. |I| displays representative samples of the obtained time series for different 
phases of an epileptic seizure. For each recording the identification of the different phases 
was done by experienced clinicians by visual inspection of the time series. 
For analysis the data of each channel were divided into consecutive segments of length 
= 4096 with a 50% overlap. The length corresponds to a duration of 24 s per epoch and 
was chosen to achieve at least quasi-stationarity for each segment according to |^ . 



To apply the method described above to these data the spectrum of each segment has to be 
estimated and for each channel a reference spectrum has to be found. 
To estimate the spectra Si the periodograms 

N 

PerM) = ^/N\J2^^J(^''^''^T (6) 

were smoothed: 

Si{uk) = '^WuPer{uJk-u) ■ (7) 

u 

As smoothing kernel w„ the Bartlett-Priestley window was chosen 

rC(l-(H)2) : \u\<b 

Wu= { n 17 (8) 

[ : \u\> b . ^ ^ 

Variance and bias of the estimator depend on the width i? = 26 + 1 of the smoothing window 
and the structure of the true spectrum. To find an appropriate value for the window width 
the spectral entropy 

H = -J2SM\nSM (9) 
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was calculated as function of B for different segments. Fig.|^a shows the plot obtained for 
one segment of an EEC The graph can be divided into two regions. For small values of B 
the fluctuations of the periodogram are suppressed insufficiently: because each summand in 
eq. (9) is a convex function of S{uJk) the spectral entropy is underestimated. In this region 
the estimated spectral entropy increases fast with increasing B. For large values of B there is 
an area of small increase where the periodogram is oversmoothed. Since information about 
the structure of the spectrum is lost, in this region the spectral entropy is overestimated. 
Fig.pD shows the corresponding plot for one realization of an AR(2) process 

X{t) = aiX{t - 1) + a2X{t - 2) + e{t) (10) 

with 



ai = 1.3, aa = -0.75 and e{t) E WAfiO, 1). 
of length N = 4096. This process describes a damped hnear oscillator driven by white noise 



35| . The functional relationship between the parameters ai,a2 and the frequency u and 



relaxation time r of the oscillator is given by: 

ai = 2cosc<;e"^/^ (11) 
a2 = . (12) 

Because this process is linear its spectral entropy can be calculated analytically: 

H = -Y,S{uj,)lnS{uk) (13) 

k 

with 

C 

|l-aie-.--a2e^2.,|2 ' (14) 

C a normalizing constant. As the plot shows the true value of spectral entropy, denoted by 
the horizontal line in Fig.pb, is reached in the area of transition from high to low increase 
of spectral entropy. Therefore, a value of i? = 33 from this region in Fig.^ was chosen for 
the analysis of the EEC 

To calculate the renormalized entropy of the EEG spectra for each channel a reference Sr{uJk) 
has to be chosen. As mentioned before this state should be of lower energy than each other 
state of the system. Because 

-'^Sj{Uk)\nSr{uJk) > - '^Sr{uJk)lnSr{uJk) (15) 

k k 

holds for every j if 

-^Sj{uJk)liaSj{LJk) > -^Sr{LJk)^^Sr{uJk) (16) 
k k 

the spectrum of lowest spectral entropy (eq. (9)) was chosen as reference. If postictal phase 
of an epileptic seizure differed from interictal phase the corresponding segment was found in 
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the beginning of the postictal phase resulting in a course of renormalized entropy shown in 
Fig.^. Otherwise the reference was found in the interictal phase. Fig.^ also shows that the 
conventional spectral entropy as given by eq. (9) does not serve as a feature characterizing 
the ictal phase. 

To test and compare the behavior of the renormalized entropy in an epileptic seizure, EEG 
data of all patients were analyzed using the method of renormalized entropy as well as simple 
features like the variance 

^Nt^EK.-^E^m)^ (17) 

j I 

or the squared euclidean distance 

b, = Y.{S,{LOu)-Sr{LOu)f (18) 
k 

of spectra. For determination of the euclidean distance the spectra were calculated in the 
same way as were done for calculation of renormalized entropy. Also the distance was 
calculated with respect to the same reference spectrum to achieve results comparable to 
these obtained by use of renormalized entropy. A representative sample is given in Fig.^. In 
Fig.^ the EEG recording, in Fig.^ the course of the estimated variance and in Fig.^ the 
course of the squared euclidean distance obtained for this recording are shown. By means of 
these simple characteristics a reliable identification of the different phases (interictal, ictal 
and postictal phase) of an epileptic seizure is not possible. The squared euclidean distance 
which was chosen as alternative and more elementary distance measure of spectra fails 
to distinguish between the different phases. The variance detects the ictal phase but miss- 
classifies a later postictal segment of the EEG. By way of contrast, the course of renormalized 
entropy reveals a temporary strong decrease only within the ictal phase. 
To investigate the spatial behavior of renormalized entropy, for each patient up to eight 
channels, corresponding to recording locations of different distance to the epileptic focus, 
were analyzed. In all patients the value of renormalized entropy within the ictal phase was 
found to decrease with decreasing distance of the recording location to the epileptic focus, as 
shown for a representative example in Fig.^. Thus, a technical device for localizing epileptic 
foci, based on the concept of renormaized entropy, is imaginable. 

IV. CONCLUSIONS 

The method of renormalized entropy, formally introduced to quantify the complexity of 
the different regimes of a dynamical system, has been applied to invasive EEG recordings of 
ten patients suffering from temporal lobe epilepsy. 

In all patients the course of renormalized entropy obtained for recording locations nearby 
the epileptic focus shows a strong decrease in the ictal phase of an epileptic seizure with 
respect to the interictal or postictal phase. Because the strength of this decrease depends 
on the distance of the recording location to the focus not only a discrimination between the 
different phases but also a localization of the focus seems to be possible. 
The method makes exclusively use of the spectral properties of the time series under consid- 
eration and therefore human interaction is restricted to the choice of the spectral estimator 



6 



to be used. 

Putting it altogether the concept of renormalized entropy seems to be a promising candidate 
for clinical applications like seizure detection or localization of epiletic foci. 
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FIGURES 



FIG. 1. Invasive EGG recording of an epileptic seizure (a): segments of the (b) interictal phase, 
(c) ictal phase and (d) postictal phase. The vertical lines in (a) denote the beginning and the end 
of the ictal phase. The measured local field potentials are shown in arbitrary units. 

FIG. 2. Estimate H of spectral entropy versus width B of the smoothing window obtained for 
a segment of an EEG (a) and an AR(2) process (b) of length N=4096. The horizontal line in (b) 
denotes the true value of spectral entropy. 

FIG. 3. Course of spectral entropy H and renormalized entropy /S.H obtained for a recording 
location in the epileptogenic area. Vertical lines denote the beginning and the end of the ictal 
phase. 

FIG. 4. Course of (b) estimated variance a^, (c) euclidean distance D and (d) renormalized 
entropy Ai7 obtained for the EEG shown in (a). Vertical lines in (a) denote the beginning and 
the end of the ictal phase. 

FIG. 5. EEG and course of renormalized entropy obtained for recording locations of different 
distance to the epileptic focus: (a) location nearby the epileptic focus, (b) location of smallest 
distance to the epileptic focus, (c) location on the contralateral hemisphere. Vertical lines denote 
the beginning and the end of the ictal phase. 
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